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Abstract 

Geometric discretisation draws analogies between discrete objects and 
operations on a complex with continuum ones on a manifold. We gen- 
eralise the theory to the cubic case and incorporate metric, by adding 
volume factors to our discrete Hodge star and then by modifying our 
inner product which leads to the same result. 



1 Introduction 

Geometric discretisation (GD) ^ [2] appears quite complicated, at first, using 
unfamiliar maps and objects but nothing could be further from the truth. Sim- 
ply put we can translate from continuum objects to discrete ones, and back 
again in such away that not only can we make the discrete ones as close as 
we want to the continuum ones, but we have discrete operators which satisfy 
the same identities as their continuum counterparts. We have d 2 = 0, Stokes' 
theorem, the Leibniz rule and more; meaning that the discrete theory mimics 
the continuum one to a remarkable extent. Here we are interested in the theory 
itself, as opposed to applications of which there are many. In particular, we 
introduce a cubic Whitney map |3] and metric, needed for the discrete Hodge 
star which is after all not purely a topological object. 

The basic structure of the theory, using the Whitney and de Rham map, to 
translate from the discrete the the continuum and back again, is the same as 
considered by Dodziuk ^ . GD uses a subdivided space though, in order to have 
a discrete Hodge Hiptmair [H] considers the discrete Hodge star using finite 
elements which map also to dual spaces and satisfy 8 = *d*, where he refers 
to cubic work done by Nedelec [H] which I have only recently come across. My 
focus has been on application to lattice field theory [7] , though the relationship 
with finite elements has also been of interest. 

GD deals with operators, as well as functions, in such a way that the identies 
and rules which they obey hold. Notably we have 

• Stokes' Theorem (Gauss' Law in electromagnetism) 
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• The Leibniz rule (The product law for differential forms) 

• A discrete Hodge star * where * 2 = I and 8 = 

• d 2 =0 

• A skew symmetric wedge 

and DO NOT have associativity of the wedge. The presence of a Hodge star 
with the associated properties is the most significant feature of the system. 

There are two aspects to our work here. First, we introduce a cubic Whitney 
map and then tackle convergence. We find that by the addition of volume 
factors to the Hodge star we can demonstrate convergence. This is not entirely 
satisfactory since the inner product is intimately related to *; we cannot alter 
one without the other. That being the case, we look at how a natural inner 
product leads to the introduction of exactly the volume factors which we had 
anticipated above. 

We have also been considering the implications of this theory to lattice field 
theory, where the lack of a discrete Hodge star has been a problem [8] but have 
been aware that related work has been ongoing from engineering E] ■ The 
possibility of fruitful cross-fertilisation is very much on the cards and something 
which we have always been interested in developing further. 

In short we begin with a review of GD before going on to our new work 
where: 

1. We develop a cubic theory. 

2. We demonstrate converge using a heuristic involving "natural volume fac- 
tors" . 

3. We show that convergence, using a new inner product which retains the 
relationship of the Hodge star to the inner product, leads to precisely the 
same factors as we used for the heuristic. 

We finish with a brief discussion of our current and future work. 

2 Geometric discretisation 

Geometric discretisation (GD) is a discretisation scheme based on a corre- 
spondence between discrete objects and operations on a complex, K, |101 ITT] 
with continuum ones on a manifold, M, 1% which captures topology 

The de Rham, A K , and Whitney map, W K , play a central role since they 
allow us to move from continuum to discrete and back again g] whilst main- 
taining topology 1 . For the moment, we note that this provides a very clean 

x Both these operators commute with d, so if / is exact or closed in the continuum, it is 
also in the discrete. 
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structure and allows us, for example, to induce a natural discrete wedge from 
the continuum one 

X A K y = A K {W K {x) A K W K (y)), 

which inherits skewsymmetry and the Leibniz rule; though not associativity. 

The other key idea is that of a subdivided space B, containing both the orig- 
inal, K, and dual, L, spaces, which allows us to introduce a discrete Hodge star 
operator £Q. This has the property that ** = 1 and S = with appropriate 
signs, which is of interest 8 . For example, we are able to capture chirality in 
the Dirac-Kahler formalism as a result [7|. 

In summary we have the following in GD, with the noticeable exception of 
associativity for the wedge: 

• An exterior derivative which maps from p-cochains to (p + l)-cochains 
with d 2 = 0. 

• A wedge product, with which we can take the product of a p-cochain and 
a g-cochain to get a (p + q)-cochain. This has the properties that 

— Skewsymmetry:x p A y q — (—l) pq y q A x p . 

- Leibniz Rule:d(x p A y q ) = dx p A y q + (-l) p x p A (dy q ). 

• The Hodge Star: This duality map associates an (n — p)-cochain in L to 
each p-cochain in K, capturing -k-k — I and d — *d-k with appropriate signs. 

For associativity we get 

(** A K y q ) A K z r = (£±f±£) * P A* (y q A K z r ). (1) 

Given a triangulation 2 of our space we can translate from differential geom- 
etry to our discrete structure. We start by looking at how forms are projected 
onto the triangulation and back before moving on to various operations. 

2.1 de Rham map 

p-cochains, C P (K), on the triangulation are the discrete analogies of p- forms, 
J7 P (M), on the manifold. The p- forms can be projected onto the p-cochains using 
the de Rham map, A K , which involves integrating them over the associated p- 
chains. In other words we have A K : fl p (M) — > C P (K) defined as 

<A K uj p ,a p >= [ cu p , 

i 

which has the property that dA K = A K d K (Stokes' Theorem). 

2 We do not worry about how to construct a triangulation of the space we are working with 
though we do know that they exist for the cases we are interested in since Rado 13 proved 
this for compact spaces. In fact we know any differentiable manifold can be triangulated !14l . 



3 




d 



P+l 



c 



A 



w 



A 



A 



W 



P 



d 



P+l 



Figure 1 : d, A K and W K provide a commutative diagram if we restrict ourselves 
to the image of the Whitney map (a.k.a. Whitney elements), a finite dimensional 
space of functions. 



If we take the standard triangle [0, 1,2], shown in Fig[21 we can see that what 
happens explicitly. Since this is in 2D we have 1, dx, dy and dxdy as our possible 
differential forms, all with possible function coefficients /. The possible cochains 
are [0], [1], [2], [0, 1], [0, 2], [1, 2] and [0, 1, 2]. 1 is a 0-form and so is mapped onto 
the vertices. This means that / is mapped to /([0]) + /([l]) + /([2]). Similarly 
fdx is mapped to (J [Q1] fdx)[0, 1] + (jj 2] fdx)[0, 2] + (J [12] fdx){l, 2]. 

2.2 Whitney map 

The Whitney map is the complimentary operation, from p-cochains to p-forms. 
In order to introduce this we need the barycentric coordinates, /V s - Given an 
n-dimensional complex we have fig, ■ ■ ■ , Mn defined on each n-simplex which have 
the property that: 

• Mi(N) = %• 

• J^. Hi(x) — 1 for all x inside the triangle. 

• fit = outside the triangle. 
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Figure 2: The standard triangle. 



4 



In our standard triangle, with vertices with coordinates (0,0), (1,0) and (0, 1), 
we can define 

Ho = 1-x-y 
A*i = x 

A*2 = y 

which satisfy the conditions necessary. 
We can then define W K as 

W K [v , ...,v p ] = ^p!(-l)V i d^ A • • ■ dfi 1 ■ ■ ■ A dm 

i 

where d[i l denotes that dm is emitted. With this we have that dW — Wd pflllSj 
and A K W K = I. This leads to a commutative diagram of sorts, Fig. ^ since 
we can starting from a form, map it to the triangulation and then act on it with 
d or act on it with d before taking our approximation; either way we get the 
same result. It is a "commutative diagram" since W K A K is not equal to the 
identity, which means that not all routes are equivalent. 

2.3 Wedge 

Having established the basic mechanic of using W K and A K , to map to and fro 
between our spaces, we can induce a discrete wedge product f\ K on the discrete 
side: 

x p A K y q = A K (W K (x p ) A W K {y q )). 

This is both distributive and anti-symmetric but not associative 3 , as can be 
seen from Eqn^ 

2.4 Hodge star 

The Hodge star is the jewel in the GD 4 crown. There are various problems with 
having a discrete star satisfying both = I and 8 = -kd* type behaviour as 
discussed by Rabin J5J . These can be resolved by working in a subdivided space 
and saying that * maps from the original triangulation to a dual space. The 
advantage of this can be seen from the Figures 13151 

First from Fig. [3]we see that with a dual space we have a trivial identification 
of original and dual objects and so capture ** = /; while in Fig. 21 we see that 
without a dual space, we do not return to where we began and so do not have 
** = /. Finally, in Fig. [5] we see that without the use of a subdivided space we 
do not get S = *d* but end up linking vertices to vertices "two units away" . 

Aside from associativity of the wedge, we lack a Whitney map on the dual 
space. Partially motivated by this, we develop a cubic version of the theory 
next, whose dual space is also cubic. 

3 It is an anti-symmetrisation like 1161 and 1171 

4 This was developed by Adams 
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Figure 3: A dual space leads to ** = I type behaviour. 




Figure 4: Without a dual space * doesn't know which vertex to map back to 
and so end up mapping to all of them. 

• • 1 • • 



1 . I 



Figure 5: Without a subdivided space we do not get = ★ 6* type behaviour. 
Initially * maps the vertex to the edge, before d maps to its edges. Finally the 
second * maps these vertices to edges which results in a line which is twice as 
long as it would have been if it initial vertex had been acted on with d. The 
use of a subdivided space halves the lengths and thereby fixes the problem. Of 
course multiply both sides by * leads to d = *5*. 
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Figure 6: Binary notation applied to a square, where vertex number expressed 
in binary leads leads to the x and y coordinates of that vertex. This makes 
it very easy to know which combinations of vertices are edges since they must 
differ by a power of two. In 3D for example we can tell that the edge [2,3] is 
connected to [6,7]. 



3 Cubic theory 

We begin by introducing our notation before looking at the operators d and d 
in this language and moving on to the de Rham and Whitney maps, which we 
show to have the desired properties. 

Note that the dual to a cube is also a cube which means that we also have 
a Whitney map from the dual space which was not the case with simplices. We 
unfortunately were unable to develop a generic method to determine Whitney 
maps, which is the tricky part to generalise, for arbitrary cells. 

3.1 Cubic Notation 

We describe the various vertices, edges and faces, taken from Fig. |HJ in our new 
notation below: 



In n-dimensions, our general object, for the unit hypercube, is thus [clq, . . . , a„], 
where a, can be 0, 1 or Xi. 



Vertex 
Vertex 1 
Vertex 2 
Vertex 3 

Edge from to 1 
Edge from to 2 
Edge from 1 to 3 
Edge from 2 to 3 
Face "[0,1,2,3]" 



[0,0] 
[0,1] 
[1,0] 

[1,1] 

[o,y] 

[x,0] 
[x,l] 

[i,y] 
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3.2 Cubic Operations 

Wc define D = d cube as follows: 

D[a ,...,a n ] = y^SiPi[ao, . . . ,a n ], (2) 

i 

where Sj = if a.j = Xi, s» = +1 if dj = 1 and Sj = — 1 if <Zj = 0. Pi equals 
(— 1)", where n is the number of a/s equal to Xj for j < i. 
Next we introduce the boundary operator: 

d : [a , ■ ■ ■ ,a n ] = ^Pi([a , • • ■ ,a n ]- [a , ■■• , 0, • • • ,a n }), 

i 

where we have replaced the ith slot with a and a 1, the second one having 
opposite the orientation, is the same as for D. In other words 

d[x] = [1] - [0] 

and 

d[x,y] = -[0,y} + [l,y] + [x,0}-[x,l] 

which is what we expect the boundary operation to be. To test it more thor- 
oughly we apply it to [x, y, z\. 

d[x,y,z] = -[0,y,z] + [l,y,z] + [x,0,z] - [x, l,z] - [x,y,0] + [x,y, 1] 
d 2 [x, y, z] = [0, 0, z] - [0, 1, z] + [0, y, 1] - [0, y, 0] + [1,1,*]- [1,0, z} + 

[1, y, 0] - [1, y, 1] + [1, 0, z] - [0, 0, z] + [x, 0, 0] - [x, 0, 1] - [1, 1, z] + [0, 1, z}+ 
[x, 1, 1] - [x, 1, 0] - [1, y, 0] + [0, y, 0] + [x, 1, 0] - [x, 0, 0] + [1, y, 1] - [0, y, 1]- 
[a:, 1,1] + [*, 0,1] 
= 0. 

The boundary map is clearly taking us to faces of the object it is acting on but 
from this we can see that the orientation of the various factors is being dealt 
with properly too. 

We want to show that d 2 = in general. We can see that d maps an m-cell 
to two oppositely oriented pieces. We say that the first d sets a,k to either or 
1 and the second the same to a/. We denote this as 

m „m— 1 „m—l 
Off = (Tq — (Tj 

in short hand, where we mean sum over i where at = by er™ -1 - Then 

oft m _m— 2 _m — 2 _m — 2 . _m — 2 
d cr =(T 00 -C7 01 -C7 1Q +C7 n . 

In fact we have two cases to consider depending on whether k < I or k > I. The 
idea is that if k < I then we get the same sign factor whether the first d removes 
k or the second. This is not the case with I since the pi factor changes. Thus we 
get two oppositely oriented versions of the same term, which cancel. The same 
argument applies if k > I. 
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dW K = W K D 



3.3 de Rham and Whitney 

We next introduce A K and W in the cubic framework, showing they satisfy 
the desired properties 

• A K W K = I 
• 

• DA K = A K d 

The Whitney map is defined as 

W K [a ,...,a n } = W K [a ]W K [a 1 }...W K [a n ], (3) 
where, for edges [0, Xj] of length hj, we have 

W K [0] = (4) 

W K [1] = (5) 

W K [ Xi ] = (6) 

From EqnO the general case consists of products of terms of the form 
W K [a t }. W K [a,} maps to = ^ = 1 where a* = x* = 0, to f*- = ^ = 1 

where a, = x, = hj and since W-^Xj] = we get 

o 

when aj = Xj. Thus A K W K — I in general. 
From Eqn[21and Eqn|3]we see that 

W*D[ao, . . . , a n ] = SiPiW K [a Q ] . . . W K [a„] 

i 

and 

dW K [a ,...,a n ] =Y. wK ^--- dwK ^--- wK ^ ( ? ) 

i 

= s lPl dW K [a t ]W K [a ] . . . W* [a,] ...W K [a n ] , (8) 

i 

Thus 

dW K = W K D. 

Finally, from Stokes' theorem we see that A K d acting on a p-form 



^[of 1 ] 
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is the same as 




= dA K w". 



For this to hold in the cubic case we need d and D to be compatible, as they 
are in the simplicial case via the incidence matrix. 

Looking at d we can see that the sign, of incidence matrix elements, is 
determined by whether we introduce a or a 1. This is also the case in D. We 
get a (-1) contribution if we introduce or remove a 0. The remaining factor is 
the pi one which also occurs in both d and D. So the incidence matrix |14| 
associated with d induces a D which is the same as the one we have that been 
using. Since D is compatible A D = dA K follows. 

Note that since A K W K — I we also know that 

D 2 = (A K W K )D 2 = d 2 (A K W K ) = 0. 

Thus we have 

1. Introduced a cubic notation. 

2. Defined appropriate boundary and coboundary operator. 

3. Developed a cubic Whitney map. 

4. Shown that these satisfy the desired properties. 

Note that if we map dx onto the standard triangle then due to the diagonal 
edge, [1,2], we introduce a dy component when mapping back using W . This 
does not happen in the cubic. 

Next we use this formalism to incorporate metric into geometric discretisa- 
tion and show that image of the discrete Hodge star converges to the continuum 
one. 

4 Incorporating metric via heuristic 

We have treated GD primarily as a topological theory but ★ has metric depen- 
dence which must be considered also. Since we can only compare two similar 
objects, we map continuum ones to the triangulation and back again, using 
W K A K , and note the difference. So when looking at functions, or forms, we 
consider 

/ - W K A K f, 

and its dependence on the size of the discrete cells used. 

In ID, for example, we can see that via the de Rham map, f(x) goes to 
/(0)[0] + /(1)[1]. W then maps this to /(0) + (/(l) - f(0))x; a piecewise linear 
approximation. For fdx we get J, Q .-, fdx on the complex which W maps to 
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(/[oi] fd x )df J >0- Since fio goes from to 1 as x goes from to L, the length of 
[01], we get 

W K A K (fdx) = w[ f<l,,, — ' 



- [01] L 

or the average value of / along the edge [01]. 

What we find is that the approximation made involves taking the average of 
the continuum object for forms and piecewise linear approximations for func- 
tions. The de Rham map takes the integral over a triangle, say, while the 
Whitney map results in dividing this by the volume so 



v 



f<Fx -» ( J f(Fx)o* - J ^ gp - /averaged. 

We can analyse * using this picture. If we start with a p-form, we integrate 
it over a p-simplex. This leads to a p-cochain in our discrete structure. We can 
act on this with *, which maps it to an (n — p)-cochain which we can be mapped 
back to continuum space with Whitncy(if we had a dual Whitney map that is, 
which we do for the cubic case). The problem then is that the Whitney map 
divides by the "wrong" volume; we no longer have an average: 

r ( f fd p x) 

f<Px -> (J ^ fJPx)o> - J £lU-r X + f^° dPx - 

The idea then is to introduce volume factors to * which fixes this. We say 
that 

Vokr"~ p 

'O 7} 

*NEW& y — > 77 : „ * cr F . 

VOlCjP 

Using this we have anticipated the problem with W so that now, when we divide 
by the "wrong" volume factor, it simply cancels with what we have, leaving the 
correct term to get the average. This also has the property that when we look 
at ** the factors cancel out, as they should. 

Once we have this, we can see that S = -kd-k and the Laplacian should also 
converge. For the first case we want W K -kd-k A K to converge to the continuum 
case. We can rewrite this as W K * A K W K d-k A K since A K W K = I, but this is 
{W K *A K )d(W K *A K ) since dW K = W K d. So if we have W K *A K converging 
to *, which we do with the addition of the volume factors, we have the desired 
result. The Laplacian follows once we have 6 converging since we can move d 
through A K and W . Note that we only have these result for the cubic case 
since we do not have * acting on the dual space otherwise. Dodziuk 0] has 
previously shown convergence of functions, forms and various operators though 
his system was did not involve a subdivided space. 



5 Incorporating metric via inner product 

The previous modification is nice, since with it we have our discrete * converging 
to the continuum one BUT it is against the spirit of GD. By adding an ad hoc 
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term we break the relationship between the inner product and *. We desire our 
discrete theory to be as close as we can to the continuum one. We look at the 
inner product and what we can do with that next. 

5.1 The inner product is the star 

Any linear map on a vector space can be expressed in terms of the dual basis 
list elements. If we have any p-form A we can wedge it to an (n — p)-form fi to 
get an n-form fda = A A [i. This is a linear map so / is uniquely determined. 
Since it is linear we can also express it as < *A, [i >= / where *A belongs to 
the dual space. This * is the hodge star and we can see how closely it is related 
to <, >. We cannot change one without changing the other. 
We had in our original formulation that 

(e i ,e?)=6 ij 

which means that the inner product of an edge with itself is 1 whilst we'd expect 
from the continuum case to get its length. This can be achieved by using a new 
inner product 

(e\e 3 ) = J W(e*) A*W{e j ) 

as used by Dodziuk 0|. We want to investigate the effect this has on *. Wc 
have introduced metric information into our system and ideally this should filter 
through the system leaving all the properties which we are happy with whilst 
sorting out the metric dependence of *. 

5.1.1 Determining * 

Our inner product is given by 

< e B , eP B >= J W B (e B ) A *W B (eP B ), (9) 

Note that e l B and e J B both need to be of the same dimension since we need an 
n-form on the RHS for it to be nonzero. We then define *k as follows: 

<* K e l K ,e%>= Jw B (e l K )AW B (eZ), (10) 

where using the fact that e K and e™ can be expressed in terms of elements of 
B. 

Now we write 

* K e l K = h l m e?. (11) 

We can using © and (|10fl determine h\ n which determines our matrix for 
the hodge star operator. 
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Firstly we define two more matrices A and B which allow us to move from 



basis elements of K and L to B respectively. So 

4 = 44, (12) 

e£ = Bfe? B . (13) 

and (|10f> can be rewritten as 

<* K e l K ,e%> = <h l e° B ,e%>, (14) 

= h l a B?BJ ><e B ,e B >, (15) 

= h l B?Bf J W(e B ) A*W(e B ). (16) 

I ij = J W(e%) MrW(eP B ). (17) 

<* K e l K ,el l > = h l B°PiBJ\ (18) 

= h l X om . (19) 

Also using l|ll)|l we have that 

<* K e l K ,el l > = J W(e l K ) AW(e%), (20) 

= Jw(A{e B )AW(B™e B ), (21) 

= A\B™ j W(e B )AW(e B ). (22) 

^ = Jw(e B )AW(e B ). (23) 

<^e^,e2 l > = X-I'^ir. (24) 

= S 1 '™. (25) 

Since l(T5 )l =l(25 )l we know that 

h l Q X om = S lm , (26) 



which is just a matrix equation. If we determine X 1 we can right multiply by 
this to get h which is *k- We can get *l using a similar calculation. 

5.2 Cubic case 

Given an inner product, the Hodge star operator can be determined, as we have 
just seen. We consider here the cubic case 5 and show that volume factors emerge 

5 We have investigated the simplicial case computationally but the nature of the volume 
factors is unclear due to its more complicated form. We resulted with ** diagonal but not 
equal to the identity; this could be normalised but the metric dependence, which was the 
objective, was not clarified. 
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Figure 7: The reference square. 



in precisely the manner we expected from the heuristic. We take a simple one 
square system which is subdivided as our reference system. 

We begin by determining the Whitney map of the various cells which we 
have. For example 

W B [0] = ^(a-x)(a-y) A 
= ^{x){a-y) B 

= M x )(y) C 

= ^(a-x)(y) D 
W B [7] =-i(a-x)y C 

W B [0123] = £dx A dy A, B, C, D 
Once we have this we can look at <,>k where 



y p K >k= J ^(b^)a# b K). 



If we let the e l be the basis list elements of K then we can define 

I ij =< e\e> > K . 

Note that 

• I is diagonal otherwise e l and are in different squares of B and so their 
images under W B ( though not under W K ) do not overlap. 

• The various integrals that appear are the same since we are always inte- 
grating either x 2 from to a or something which can be expressed as this 
with a change of variables. 

• The dx and dy integrals are independent and both of the form mentioned 
above. 
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The are only two possible cases which lead to the same result: 



/ 

Jo 



a 3 



x 2 dx = [-]g 

o 

° t 

3 



f a f a 

/ {a-xfdx = I y 2 {-dy) 

J0 J a 



So any of the integrals which we have are equal to 

In the case of vertices we have two such integrals leading to factors 
with a symmetry term F to specify how often a term occurs; this is 4 for vertices 
since each vertex occurs in 4 squares, 2 for edges and 1 for faces. For edges we 
only have one integral, since we don't get x 2 y 2 terms but x 2 dy ones instead, 
and for faces we get no such terms since the integral is simply dxdy. B also 
introduces a factor, as in the simplicial case; though here we get 2 P instead of pi 
since the volumes involved are different 6 . Finally W K has a term associated 
with it so for vertices we get 

1 x 4« 6 
<[0],[0}>k- * ' " 



for edges 



and for faces 



< [01], [01] > K = 
< [0123], [0123] > K - 



9a 4 ' 

4 x 2a 4 
3a 4 

16 x la 2 



For L we similarly get 

< [0123], [0123] > L = 
< [01], [01] >l = 

<[6],[6]>i = 



1 x 4a 2 
4x2 



3 

16 x 1 



We only mention the results for particular vertices, edges and faces but they are 
all the same 7 . 

To determine *k we use 

<*KX p K ,y n L - p > L = J W B {BxV K )NW B {Byl-V). 

We have the LHS except *k so we need the RHS next. Note that since both 
<, >k and <, >l are diagonal, we only get contributions when y = x. 

6 This is related to maintaining A K W K = I. 

7 The F factors are all the same since we are interested in complexes without boundary in 
which case every vertex has F = 4, every edge has F = 2 and every face has F = 1. 
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We have two B factors on the right hand side which leads to a factor of 4. 
If we are dealing with a vertex then we wedge it with a face. B of the vertex 
gives 1, while B of the face gives 4. For edges we get 2x2 since we wedge two 
edges together and for a face we wedge it with a vertex leading to a 4 x 1 = 4 
factor. We look at the various cases next. 

Performing the RHS integral for the vertices, if we have x = [0] then y = [0] , 
we get: 

= 4 f — r(a — x)(a — y)dx A dy 

J A a 

a 4 2 2 
= 1. 



For the edges we get: 



= 4 J W B [01] AW B [03] 



1 a 2 a 2 



A 



And for the the faces: 



/ 



W B [0123] A W B [0] 



I A a 4 2 2 



= 1. 



The dual calculations are the same apart from sign factors which arise when 
the W K (x) and W K (y) terms are flipped. So for the edges case we get a (— 1) 
sign. 

Determining -kx is now trivial since both sides of the equation are diago- 
nal. The ith diagonal element of is simply the ith diagonal element of the 
RHS(which are all just ± 1) divided by the zth diagonal element of <, >l which 
only depends on whether dealing with vertices, edges or faces. So on 0-cochains 
have 

r , ft 2 r", 

**[0] = jg[0] 

on edges we have 

**[oi] = |[rii] 

and for faces we have 

*k[0123] = ^[0123]. 
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We can get *l in the same way: 

*l[6] 

*l[01] 
*l[0123] 

With this we get 

which is diagonal and proportional to the identity. 

We can then introduce normalisation factors on the RHS — as is done in the 
simplicial theory — by simply adding a factor of | to the definition of 

< * K x p K , V r P >l= 1 1 W B {B X V K ) A W B (By n L - p )- (27) 

Using this we have *k*l = H as required. 

Here whilst we get = I type behaviour, it can be improved. *k acting on 

a vertex introduces a factor of i, after the above normalisation has been made. 

6 1 

When *£, acts on this it multiplies it by 6 so we get the desired ** = I. Since 
does not occur in the continuum we should normalise at this level instead which 
can be done by using appropriate p dependent factors in Eg 1271 

We can now determine convergence results since we can see how the length 
scale a appears in the various operators. In short they appear as you would 
expect them to using the heuristic which means that that S and the Laplacian 
converge to the continuum case. 

Firstly we look explicitly at what happens for S = -kd-k in 2D. In the contin- 
uum case we know that: 

1 . -kd-k f = -kdfdxdy = 

2. -kd -k fdx — -kd(fdy) = *^dx A dy = |^ 

3. *d*fdxdy = kdf = *(%dx + ^dy) = ^dy - &dx. 

We can now look at what happens in the cubic case and compare W K -ki J di J ^K^ K 
with the above. 

Note that d = -kd-k, with some sign factor which means that -kd-k applied to 
[xy] gives us its boundary, upto sign. For example, if we take the edge [0, y] we 
map —fi onto it, where fa is fdxdy integrated over the face i. The adjacent 
face will also contribute to this edge with a + fj term. Face i is a away in the 
x direction from face j and so we effectively get an f(x + a) — f(x) term which 
after Whitney will give 

f(x + a)-f(x) 
a 



- -> 

= ^[0123] 

= ±i-J 
64 
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f(x) 


f(x+a) 


• 


• 



Figure 8: Oppositely oriented adjacent regions lead to f(x + a) — f(x) terms. 



In this way we can see how we agree with the continuum result. The 1-form case 
follows similarly with the boundary now vertices and the 0-form case is trivial. 

As we mentioned before, the volume factors which arise are precisely those 
which we required from our heuristic ( Sec^J) which can be generalised to higher 
dimensions. First we note that this is the case in 2D since * applied to vertices 
leads to a 2 factor which is the volume of the cell mapped to divided by that of 
the cell acted on. For edges we get a factor independent of a which again agrees 
with this picture as does the factor for faces. 

In general we have a 3D ~ 2p /a 2D for p-cochains. In two dimensions this means 
that we have a 6 ~ 2p /a 4 so a 2 for vertices, 1 for edges and a~ 2 for faces. Similarly 
in three dimensions we have a 9 ~ 2p /a 6 or a 3 for vertices, a for edges, a -1 for faces 
and a -3 for cubes. In general the heuristic expression is the volume of the space 
mapper to, a D ~ p , divided by the volume of the object being acted on, a p . Since 

3D-2p D-p 
^ _ n D-2p _ ^ 

a 2D a p ' 

the two agree. 

We had already seen that the heuristic suggested a modification that could 
be made to * so that it, along with the coderivative and Laplacian, would 
converge to the continuum result. This problem with this was that it broke the 
relationship between the inner product and *. In order to retain this we made a 
natural modification to the discrete inner product to see what effect this had on 
-k K and found that this resulted in adding precisely the factors which we wanted 
from the heuristic. As a result, we can use the various convergence results which 
we had before, except that now the relationship of * A and the inner product is 
preserved. 



6 Conclusions 

We have thus introduced a cubic Whitney map which we have used to demon- 
strate convergence of our discrete functions and operators to their continuum 
counterparts. This involved the introduction of some modifications for the 
Hodge star, either via a heuristic, which was unsatisfactory due to the rela- 
tionship of the inner product and *, or via a new inner product (the one in fact 
used by Dodziuk). 
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With metric and convergence results, applications become a possibility. The 
work, for example, of Bossavit Hiptmair |S] and Nedelec [H] are naturally 
of interest though within a different context, namely finite element. GD also 
involves complexes, with finite dimensional subpaces within them, but does not 
use variational arguments for existence theorems for example. 

Our current focus is on application to lattice field theory where following 
Rabin [5] we have shown that within the Dirac-Kahler formalism [T^], where 
Nielsen- Ninomya is n °t applicable as shown by Becher 120], we have chirality 
(* in this picture) whilst avoiding degeneracy. 

Work has also been done extending the mathematical structures which GD 
deals and should soon follow. 
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